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16 Abstract 

A two-dimensional hybrid stress finite element analysis is described which has been 
used to study the local stress field around delamination cracks in composite materials. 

The analysis employs a crack tip singularity element which is embedded in a matrix inter- 
layer between plies of the laminate. 1 ults are given for a unidirectional graphite epoxy 
laminate containing a delaminutlon en .ting from a surface crack through the outside ply. 
The results illustrate several aspects of delamination cracks: (1) the localization of the 
singular stress domain within the interlayer, (2) the local concentration of stress in the 
ply adjacent to the truck, (3) the nature of the transverse normal and interlaminar shear 


stress Mstributions. and (4' the relative magnitudes nf Kj and Kjj associated with the 
delamination. A simple example of the une of tht analysis in predicting del’ mination crack 
growth is demonstrated for a glass/epoxy 1: ntnate. The comparisons with experimental 
data show good agreement. 
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"OR I. WORD 


This report describes a portion of the analytical and experimental 
work for the period 1 May, 1976 through 30 April, 1977. The study was 
conducted for the NASA-hewis Research Center hv members of the 
Department oi Materia. t Sicenc id Engineering at the Mnssa< ip'^tts 
Instltuti of Technology. The Principal Investigator 'A the ovt ^il 
program is Prof. F.J. McGarry; other major art icipant. are S.5. Wanp, 
and J.F Mandell. T ho NASA-LcRC Project Manager is Dr. C.C. Chamis, 

Mall Stop 49-3. 

Efforts in this program are primarily directed towards tha 
development of f in i te el emen t anal yses for the studv of flaw growth 
and fracture of fiber composites. The wo.k described In this report 
is meant to complement the development of a three-dimensional 
analysis capability, as well as to advance general understanding of 
the subject. 
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ABSTRACT 


A two-dlmfnMlon.il hybrid stress finite elment analysis Is described 
which has been used to study the local stress field around delamination 
cracks In composite materials. The analysis employs a special cr. k tip 
element which is embedded In a matrix interlayer between plies of the laminat 
Results are given for a unidirectional graphite/epoxy laminate containing 
a dclamlnatlon emanating from a surface crack through the outside ply. 

The results Illustrate several aspects of delamination cracks: (1) the 
localization of the singular stresc domain within tin Interlayer, (2) the 
local concent rat ion of stress in the ply adjacent to the crack, (3) the 
nature of the transverse normal and interlaminar shear stress dlstributio s, 
ind (4) the relative magnitudes of Kj and associated with the delaminati 
A simple example of the use of the analysis in predicting de laminat Ion crack 
growth is demonstrated for a glass/epoxy laminate. The comparison-’ *ith 
experimental data show good agreement. 
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INTRODUCTION 


This paper i» concerned with the analysis of dc-lamln.it Ion cracks 
propagating between the plies of fiber reinforced composite laminate! . This 
general clans o problems has been of concern throughout the history of 
laminated materials for the obvious reason that the bonding ' etween plies 
is commonly the element of the material system which is the softest, 
weakest and most prone to environmental attack. In the case of f iber 
reinforced plastics, to which this paper Is specifically addressed, de- 
lamination problems have been of concern from some 01 the first appl icat ions 
involving the unwinding of missile casings and pipes, to current difficult- 
ies with del amin it ion-induced failures In aircraft structures. 

Analyses of delamination-type pt oleins have ta':en the form of very 
approximate predictions of strengtli (1,2) and damage extension rateB [3] 
or else more rigorous solutions for very simplified r.odels 1 4 . 5 ) . A full 
understanding of such problems requ'rc. an accurate solution for the local 
st re s distribution in the interply region which takes into a cou* ' the 
true geometric and material variables and structural parameters. ue 
analysis described in th 4 s paper ir. intended to provide an accurate numerical 
solution for a more ri-presentat ive model than those previously used. 

MECHANICAL MODEL AND ASSUM 1 TIONS 


The geometry of the problem involves a crack situated between plies 
which typically have fibers oriented in different directions. Each ply is 
actually heterogeneous, consisting of fiber end matrix regions, and there 
may be matrix-rich region betwen plies. Figu. I shows a micrograph of n 


_ l 




griphlte/epoxv composite containing interply matrix regions which are 
particularly prominent; other composites may have less identifiable 
Interply regions depending upon the materials and fabrication methods. 

The mechanical model used In the present analysis assumes that the In- 
dividual plies are homogeneous and anisotropic, and Includes isotropic 
Interply matrix regions, .interlayers) between each ply. The Interlayers 
may vary in thickness, but are generally on the order of 1/10 to 1/20 of 
the ply thickness, similar to Figure 1. The crack tip Is embedded in the 
matrix Interlayer as Indicated In Figure 2. Linear elastic behavio (s 
assumed for the plies and for the Interlayers. 

The crack geometries considered In the present study are B'milar to 
that in Figure 2. A surfs crack is assumed to penetrate one or more 
plies, initiating the delait. .nat ion crack. The remainder of the compsoslte 
is assumed to be well bonded and defect-free. 
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METHOD OF ANALYSIS 


Ceuer.i ! Furmu lit Ion 


The hybrid stress approach of the finite element method, pioneered 
by Pain (6,7 J, is characterized by the use of an assumed stress field in 
the element and an assumed displacement field along the element boundaries. 
The formulation of the analysis is based on the minimum complementary energy 
principle (a modified complementary energy functional 1b ui..u n the for- 
mulation of the crack tip superelement). Additional features of the method 
are r lexlbillcy of formulation, selection of elements and expedient 
achievement of Intereiement compatibility. More accurate solutions and 
faster convergence rates than those of conventional displacement models 
can be obtained. The singular behavior of the crack tip region, which 
's critical to the fracture problem, can be exactly modeled in the 
fromulaticn without an increase in the nu..iber of elements. The complex 
gem-trie variables and multiphase materials effects also are convenient- 
ly taken into account. 

The hybrid stress finite lement method procedure is specialized 
for application to the current plane crack problem by the introduction of 
a crack t lj> ? verel' : witlilr which the singular stress behavior 
is considered b ope-ly selected Lt jss fu- -ions, i" is assumed 
l less hybrid model for the problem was first introduced by Plan, 

(6], and later refined by Tone et al [8] by the us ’ -f the complex 
variable formulation of the Muskhel ishvi 1 i si r.".s f unct ion . The general 
formulation of the procedure is given here for the case of the plane crack 
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problem In a compogl t« laminate v/lth the superel ement embeddc 1 in the 
interlayer. The formulations for the crack tip super-element and Its 
surrounding regular (non-singular) hybrid elements have been described 
In detail elsewhere (6,8), and only a brief outline is given here. 

Consider a crack system of the type shown in Figure 2. The complemen- 
tary energy functional of the whole domain of the specimen (after division 
Into a finite number of discrete elements) mav be expressed as 

% * ft* I 'if 1 (1 

>- r»i 4 r 


where T ’ m la refered to the crack tip superelement and * ^\o the l 1 *' 

regular surrounding element in the given domain. Applying the variational 

energy principle to the functional 1T , one obtains 

c 



where 





represents the general form of the element stiffri-— matrix, and is cal- 
culated in different ways in the crack tip superelement and the surround- 
ing regular elements as decribed in the next section. The equivalent nodal 
force Q is defined as 
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Tht* matrices, Q , It and £ are defined In following; sections In detail. 
~n n n 

Assembling all the element stiffness matrices together, a act of 
linear equations of the form 


< % ” £ ,4) 

Is established. The dlsplacemnt field can be solved by a standard Causs- 
Cholosky elimination scheme 


% = K'' Q (» 

-vf A* AJ 

The stress parameters £ for the n 1 * 1 element arc then calculated from 

T1 

& &■ w 

The st8Goclat , 'd stres~ field a‘- the "ocaMcn (x , ) vC ."■st can 1 • 

P P 

obtained from 

S.fe. J- ”> 

where b Is a Boolean transformation. 

- n 


Crack- Tip Element F or mu 1 at Ion 

The conventional displacement model and non-singular hybrid stress 
element have difficulty handling the crack problems even In a monolithic 
material, since the use of high order polynomials as interpolation 
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function* do** not Improve the rate of convergence for tMs kind of 
problem. The reason for this 1 m that the convergence rate of the 
finite element method 1* controlled by the nature of the solution near 
the Hlngular region (8]. However, u*e of the complex variable technique 
In the hybrid element formulation permit* proper considerr . Ion of the 
stren* singularity and of higher order effect* In th crack tip region, 
and it lead* to highly accurate results with a relatively coarse mesh. 


The modified complementary energ. functioril, v , Is used for 

in 

the crack tip element. Consider a plane elasticity problem with pre- 
scribed boundary traction Tj over tb*- boundary a and prescribed displacement 

u. over the boundary s . The functional is defined in the form 
1 u 

K* j s - fu; T. J, . * $jc?, s 1-1 


The Euler equations for this functional are 





[9] 



(101 


Following Muskhel ishvl 1 i ’ s formulat ion (9), the stress and 
displacement fields in the plane elasticity problim can be expressed 
in terms of two stress functions <t (!A) arj f (3) of the complex vari le 

8 as 

% #*>>] 

-O',, *"<■*) + w '(*>1 
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and 


2 G, - qij(i)-itS ti) • V>(i) 


where both <t> < it) and i|»(i4) are analytical In the Z-pl.ine, an»‘ G • F./211-* v) 
w 1th 


»=( 3-y)/0* y) for plane stress (12) 

or ^ j ~ 4 ^ for piane strain (13) 

In order to choose proper stresses and displacements for the crack 
elen.en' - which would account for possible singularities of all otde. ( 
as w. 1) higher order terms, the following mapping function la lntro- 
di icon 


2 » W ($) ~ J l 


(14) 


with 

- 77/2 < ary £ *77/ z 


Thus, on the C-plane the stress function; <f ( C ) and , are 

analytical functions of S . Using this mapping function, Eq . (11) becomes 


o^-cr., ♦ 2 a cr, 3 * 2 { (%)/vi ’(s\i\v(')}/^ (?) 
zc '>4 -■ ■ *<\) 


By Imposing ,'e ( ctlon free boundary condition (T *T “0) on the crack 

x y 

a^ 

J. _f (T, ♦.Tj ) ds - or /<» ♦ 2. jF(*) + V (i) 


(16) 
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4< ( t. ) can b*> calculated from $( £ ) 


J/(0 * 


In constructing the super .nent stiffness matrix, $ ( f. ) 1* assus»«d 
to hr.ve the ' rm 


v l * r 

4 - 


(18) 


Thun, from Eq . (17) we have 


v (?) * - £cfej(-/jV j k^-j { 1 


»liw! e 


(19) 


b i 5 ft ♦* Am 


(non-svnv- 'tr lo rn e) 


*»i ; A 


(symmetric case) 


with the P's belie »ea* onstantb. 

Substituting liqs . (18) and (19) back Into Eqn. (15), the stresses 

and dlspl«.c.eraents ca l be explicitly expressed In terms of r as 

o-„-a0^ - %V* a c**;- 

♦ &>- (i? 1 A,,,- 

o^.o^ l • *> 

+ -1 (-n 5 ?’ ‘j4L, 

5 ^ 


(20 a.b, 
and c) 
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1 1 ... 
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From Eq i . 20 (a, 
tnd lnlcrlor displaces) 


and c; one can c:p thr boundary tractions 

as 



u - 

(V 



( 21 ; 


there f c 1. * column vector with it* component beins Pj, Bj. 

%» 

the boundary dlsplacei.nnt u shall be ass. ned in _erms of Renerallzed nodal 

• •># 

displacement q as 



A p-sbstitution of F s . (21) and (22) ino Eq . (P, and Variation of 

the func' ional n ith respect to 8 yi i.Js the crack elen.e it 
m ~c 

stlfft.jsp matrix k 

~ c 


K 


« 



(23) 


..d 


whore 


f. - 


-i 

, • / 1 

• * o 

r+ +*• Jm 


(24) 


n*i i (« T ? * T y)</ 


and S r } f ft’ L J S 

't t ^ /y, • 
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After obtaining the d lap It cement field of the syntem oy solving 

i . «ust nbled ,.obal stiffness mntr x K , he stress field in the super- 

oleuenr can be call ' it ' from Eqna. 20 ( »,b and c). The i ' ru; iu- 

* c >ltv factors K t and K WT can be relnt*' 1 to C by 

II c 



( 2 ') 


K< 'I Ivbr Id tress 1 I emei t mi mu ! it inn 


The compli lentary energy functional bo be varied for the I'ey; 1 ar 
element is giv... by 



| £ r TJ« 


whe* 


th. coiup] lance natrl/ 



-xpresseil s 


*«>.. I?‘ / 

Sii 


tVi 

for genera 1 i ' anisotropic plj of the composite. 


( 26 ) 


Along each boundary o ; he element, an assumed displacement field 
is selacten, and expressec ' .t te.ms of the nodal di..p,i.aie jnts q 



wb ! Is the Interpolation function. The stresses in the ‘nterior of the 
i lement a j expressed by undetermined stress parameters 8 
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O' s P Cx *J )/? 

' A. 

where P is choBcn to satisfy he h< ao^anoua equilibrium v.quat an 


The surface tractions, which are reltn-'d to th stre ;s components 
by T o n , can be written in th form 



t • h a 


Substitution P, ? i?7), (28), and (’0) into (263 , the functional 



where 

H * JjL R S P and * 

Taking the variation of the functional 7T i / m ^ with respect to the 
stress parameters 6 to minimize the complementary energy, the element 
stiffness matrix can Le obtained as 





V A J 


( 28 ) 


(29; 


(30) 


(31) 


( 32 ) 
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tc c u racy o f t he Ana l yai s 

The accuracy and convergence of the analysis are complicated by 
mveral unusual features of the problem and of the crack-tip superelement . 

Ah mentioned earlier, the finite element mesh must accommodate both the 
small dimonsic* the interlayer r'nickneis anti th» larger dimension! of the 
remainder of t”. specimen, a ditference of three ordern of magnitude. It 
In essential to model the Interlayer and the crack tip region accurately 
with a number of ” lemon t < across the lnterlaver tl kiess to that the high 
stress gradient within *-he inte layer may be discerned. This georeiric 
characters! Ic combines with the extreme dlffer<».ic^ in Interlayer and ply 
elastic modulus to cau®e s i > ilf leant tiMioerlca. round-off errors, so doul le 
precision inode was required for accurate so utions. 

Optimization of the mesh discretization Is also complicated by the 
differences l.i>een the crack tip superelement and the surrounding elements. 
The u' are lemei.t ,’tves an exact stress distr-h ition , and It is advantageous 
to ute c* n r cracu tip superelement as possible. On '.be other hand, 
the regular elements f. rrounding the crack tip must be sufficiently .'mail 
to giv. iccurate results in the domain beyond the superelement, which 
also is of irterest. As a re3ult, e compromise must be reached which yields 
an accurate solution both verv close to the crack tip and in the lurrounding 
region, which also gives the minimum band width for the stiffness matrix. 
Arrangement o. he mesh must alsc 1 satisfy the geometric contralut that 
the crack tip superelement be embedded in the interlayer, and that the 
number of degrees of freedom of the entire system be minimum in terms of 
computer run time. 
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A study of the accuracy and convergence of the A.ialynls has been made 
using test casts for wlilch Independent solutions are available in the 
literature. The tost cases studied Vfei e double-edge-notched plate and 
the double cantilever beam type crack geometry. Results described else- 
where (10) Indicate excellent agreement with the existing solutions for 
thi tost cases. These results, as well as convergence studies for adhesive 
crack problems similar to the present case, Indicate an accuracy within 
approximately one percent of the converged solution of Kj and for 
mesh arrangements of the type indicate! in R'gur ?. 
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RESULTS AND DI‘ JSSTON 


Str ss Distribution ' >r Unidlrect ional .uminare 


Th s *»ctton give 1 , the stress distribution for a four-ply nni- 
dlnctloi.il laminate with the fibers parallel to the* applied .‘.’•ess. The 
elastic constants for the plies given In Table 1 are tvplcal of graphite/ 
epoxv. The Interlayets between pller. have a thickness of one-tenth of 
the ply t ‘ckness and have elastic constants, also given in Table 1, 
typical of epoxy. The surface crack in this case penetrates thiough only 
one ply on the outside, and the delamina'rion cra< k has a length of three 
times the ply thickness. 

Figure 3 gives an isostx ?ss contour plot cf the distribution of the 
longitudinal stress, 0 , for * section of the laminate 1 und the de- 

lamination crack tip. The isostress contours are plotted by the computer 
and are only approximate, but give an overall view of he nature of t’lt 
stress field, lav t »! .vatu>' s of the rtresa field are evident: 

1. There is a va y localized disturbance in the stress field within 
the interlayer at the crack tip. 

?. There is a local stress concentration in the uncracked adjacent ply 
just above the delamimation crack tip. 

3. The cut ply shows some local compression alo the crack flank, and 
a gradual build up of stress beyond th del.imin ition crack as stress is 
transferred by shearing of the interlayer. 
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4. The far-fiell stresses gradually decrease through the thickness 

reflecting the bending of the specimen resulting from the nonsvmmet ry . 

The stresses in the interlayer very close to the crack tip are given 
more accurately on the log-loy plot In Figure 4. At expected from previous 
work on adhesive* [10]. The stresses very close to the crock tip follow a 
Mi7 si n'i'ilarlty Indicated by a slope of -1/2. The associated stress in- 
tensity factors, Kj and Kjj* are .0290^ an< ^ .045'' respect ivel y. The singular 
stress domain is completely • Med within the inrerlj.er. 

Figure 5 gives linear plots of the variation of n through the 
thickness at various distances ahead of the crack tip. As before, the 
severe redistribution of local stress in the ply adjacent to the delamination 
is vident. It should be noted in Figures 3 and 5 that the domain where 

0 is very high encompasses only a small part of the ply thickness ad- 

oc 

Jacent to the delamination, dropping signit leant ly within several fiber 
diameter'-. This suggests that the assumption of homogeneity within the ply 
may be questionable in this domain. This local stress concentration also 
preculudes the application of laminate theory 

Figures 6 and 7 present isostress -ontour and linear p.ots for the 
transverse normal stress, O ^ . The distribution of ^ is centered about 
the delamination crack tip and is continuous across the ply bou. darics. The 
Interlaminar shear stress, 0 given in Figures 8 and 9^18 moie dis- 
torted along tho fiber direction. 0 roaches higher values 
than does 0 for this case. 
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Crac k Growth Predicti o n: Unidir ecti onal and Croaapl le d Laminates 


Th* useufulness of the analysis in predicting crack extension uepends 
a large extent upon the choice of appropriate crack growth criteria. 

While no attempt has yet been made to select the best criterion for the 
mixed-mode crack extension, one » the most elementary criteria, the r* leal 
total strain energy release rate for composites, has been employed as an 
example. The critical total strain energy relee*- • rote criterion | 1 1 ) ~tates 
that crack extension will occur whan 

+ c j- ’ <?• <m) 


This criterion may he expressed in terms of K «.nd K as 

X a V 




(34) 


This criterion was applied by r ' ta<r 'ng valu f Kj and K fo«- 
several cracl lengths from the computer analysis. The value of applied 
stress necessary for c*ack X tension was then calculated at each delamination 
crack length from Eq . (34) and a prediction fer the delamination length,^, 

\ , the applied stress, o » was thus obtained. The composite system chosen 

OO 

for this calculation was Type 1003 glass/epoxy (3MCa.) in the plv configuration 
^0/oO/yO, due t3 the availability of experimental data for the delamlnation 
crack growth [12] as well as data for 0^ and for ar. almost Identical 

material. Type 1002 g’.ss/epo:’ ' 1], The geometry tested was -Imllar to that 
in Figure 2 except that two plies were surface cracked, one 90° and one 0®, 
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giving a surface crack depth of one-half of the total thickness. Sine*. 

the crack growth occurs In the matrix lnterlaver, the values >f anJ Kjj t , 

Were calculated from C w and i * , , by the following equations, using the 

Lc lie 

matrix elastic constants. 



I 

* 0 -*') 





The Justification for use of these relat lonrhlps is discussed elsewhere (1C). 
Table 1-c lists the values of the various parameters, and Table 1-c gives the 
elastic moduli for gla s/^povv used in the analvsls. 

Figure 10 Indicates good agreement inween calculated and measured 
delamination crack growth r-»tes under In Teasing applied stress. The 
delamlnatlon spreads in a stable manner Initially due to the decreasing values 
of Kj and Kjj with increasing crack length. The applie stresn must be 
increased to promote further crack extension until the delamination becomes 
long enough such that it propagates under constant stress. The delamination 
length in Fig. 10 is normalized by w, the experimental specimen width, 

hich was 100 t . Since the crack extension is stabl^, the calculation based on 
o 

successive crack iiitiations from different deamination lengths appears to be 
adequate, without the need for predicting the cate or magnitude of extensions at 
at each step. The calculated stress values at discreet crack lengths are simply 
fit by a curve to provide a continuous prediction. 
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SUMMARY AND CONSLUSIONS 


A two-dimenaioanl finite element analyst. 1 ? has been developed which 
gives an nccuratt solution for the local stress distribution around delaml- 
nat ion cracks in laminates. The laminate is modeled as consisting of homo- 
geneous, anisotropic plies separated by an isotropic matrix interlayer. 

A special superelement is em* dded In the interlayer at the crack tip. 

For the cases studied, the geometry consisted of a strip of material con- 
taining a surface crack intersecting a drlamlnat ion crack. 

S. ..iple results for ^ uuidl -ct tonal graph 4 te-epoxy laminate illustrate 
the fundamental aspects of the stress field, including: (1) the localization 
of he singular stress domain in the interlayer, (2) V very local con- 
centration of stress in the portion of the first continuous plv adiecent to 
the delamlnat ion, (3) the nature of the tr .nsverse normal and intet ...ilnar 
shear stress dlstrlbut losn, and (4) the relative magnitudes oi K ^ and 

associated with the delamination. An example of the prediction of delamination 
crack growth under increasing applied stress shows good agreement v.th ex- 
perimental data. 
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Table 1. Ply Plant Ic Conatanta ua**<l In Analyala 


a. GRAPH ITF/ EPOXY 


'LL 


'Tf 


'TL 


'TL 


138. Ci a i-J.OO x 10* pm; 

P._ 2 - 14.5 GPi (2.10 x 10* pal) 

5.87 CPa (0.85 x 10 ft pal) 

V • V * 0.21 
LZ TZ 


h. FIBERGLASS /EPOXY 

E - 34.5 GPa (*00 x 10 6 pai) 

4 

Ejt - E zz - 10.4 GPa (1.50 x 10* pal) 
G - 5.18 GPa (f 75 x 10 6 pal) 

kf 4 

V TL ’ V LZ * V TZ “ °' 25 

c. EPOXY INThKLAYSR 

E - 3.45 GPa (0.50 x 10* pal) 

G • 1.28 „ra (0.185 x 10* pal) 
v - 0.35 




% 
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FIGURE 1. 

MICROGRAPH OF CROSS - .SECTION OF 
GRAPHITE EPOXY COMPOSITE SHOWING 
MATRIX INTERLAYERS BETWEEN PLIES. 
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PLY 

INTERLAYER 



FIGURE 2. 

LAMINATION CRACK GEOMETRY AND F.YITE ELEMENT MLSH 
EAR DE LAMINATION CRACK TIP. 




DISTANCE THROUGH THICKNESS 


I 
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DI^TANC.. AHEAO OF CR^CK T P, x/u 

FIGURE 3. 


LONGTTUDINAL STF7.C3 CONTOURS (<r„ /a^ ) NEAR DEl.AMINA- 
TION CRACK TIP IN /Q/O/O GRAPHITE/ EPOXY. 


23 
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FIGURE 4. 

STRESSES AHEAD OF DELAMINATION CRACK TIP 
IN 0/0/0/0 GRAPHITE/ EPOXY COMPOSITE. 



DISTANCE AHEAD OF CPACK 
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FIGURE 5. 

LONGITUDINAL STRESS THROUGH LAMINATE THICKNESS AHEAD 
OF DELAf INATION CRACK IN 0/0/0/ 0 GRAPHITE EPOXY. 


DISTANCE THROUGH ThlCKNS.SS 


I 
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I 

, + O.Sl, 

ft 

1 I 



0 0.02 0.04 O.OG 

DISTANCE AHEAD OF CRACK TIP, x/l. 


FIGURE 6. 

TRANSVERSE NORMAL STRESS CONTOURS (o^cr^) NEAP 
DELAMINATION CRACK TIP IN 0/0/ 0,0 GRAPHITE, EPOXY. 
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DISTANCE AHEAD OF CRACK 

x/t, • 0.032 

x/t, » 0.75 

x/t, * 3.50 
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FIGURE 

TRANSVERSE NORMAL STRESS THROUGH LAMINATE T iCKNESS 
AHEAD OF DELAMINATION CRACK I . /0,0/0 GRAPHITE EPOXY. 



i/STANCE THROUGH THICKNESS 



0 0.02 0.04 0.06 

DISTANCE AHEAO OF CRACK TIP, x/L 


FIGL’RE 8. 

INTERLAMINAR SHEAR STRESS CONTOURS (<r ti /<T n ) NEAR 
DELAMINATION CRACK TIP IN 0,0, 0/0 GRAPHITE, FPOXV. 
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DISTANCE AHEAD OF CRAC 
x / r, « 0.032 
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FIGURE 9. 

INTERLAMINAR SHEAR STRESS THROUGH LAMINATE THICKNESS 
AHEAD OF DELAMINATION CRACK IN 0,0,0 0 GRAPHITE EPOXl. 


DELAMINATION LENGTL 



or/UTS 


FIGURE 10. 

DELAMINATION CRACK GROWTH vs. NOMINAL STRESS. 
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